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ABSTRACT 

We study the X-ray spectra of 85 BL Lacertae objects using the hardness ratios as 
given in the WGA catalogue of ROSAT sources. This sample includes all WGA BL 
Lacs with high-quality data and comprises about 50 per cent of presently known BL 
Lacs. We find that BL Lacs have energy power-law spectral indices between and 
3 with a mean value ^ 1.4. Significant differences, however, are present between 
high-energy cutoff BL Lacs (HBLs), normally selected in the X-ray band, and low- 
energy cutoff BL Lacs (LBLs), generally found in radio surveys. HBLs have steeper 
X-ray spectral slopes (a^ ~ 1.5) well correlated with aox, and anti-correlated with the 
X-ray-to-radio flux ratio and cutoff frequency, with convex overall broad-band spectra. 
LBLs, on the other hand, have flatter X-ray spectra {a^ ~ 1.1) and concave optical- 
X-ray continuum. We interpret these results in terms of different mechanisms being 
responsible for the X-ray emission in the two classes, namely synchrotron and inverse 
Compton for HBLs and LBLs respectively. The observed differences are consistent with 
the hypothesis that HBLs and LBLs are powered by essentially the same non-thermal 
engines differing mainly in their synchrotron cutoff energy. 
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1 INTRODUCTION 

BL Lacertae objects constitute a relatively rare class of ac- 
tive galactic nuclei (AGN), distinguished by their almost 
complete lack of emission lines. At present, only about two 
hundred BL Lacs are known (Padovani & Giommi 1995b) 
as compared, for example, to the thousands of quasars listed 
in the most recent AGN catalogues (Veron-Cetty & Veron 
1993; Hewitt & Burbidge 1993). 

Despite their small numbers and other somewhat ex- 
treme and non-typical properties (radio core-dominance, su- 
perluminal motion, rapid variability, high polarization: e.g., 
KoUgaard 1994; Urry & Padovani 1995) BL Lacs are nev- 
ertheless an important AGN subgroup, as the broad-band 
emission in these objects is most likely completely domi- 
nated by non-thermal processes. Their multifrequency spec- 
tra are in fact fairly smooth over a wide frequency range, 
suggestive of a common emission process, and lack the 
thermal features characteristic of the other, more common, 
AGN, like dust emission in the infrared and (possibly) disk 
emission in the ultraviolet band (e.g., Bregman 1990). 

Historically, BL Lacertae objects have been divided into 
RBLs (radio-selected BL Lacs) and XBLs (X-ray-selected 
BL Lacs), depending on the discovery band. The two classes 
showed in fact somewhat different properties, RBLs being 
more extreme, for example, in their optical polarization, 
core-dominance, and variability properties (e.g., Jannuzi, 
Smith & Elston 1994; Perlman & Stocke 1993). This division 



was clearly not satisfactory, as it was not based on intrin- 
sic physical properties but solely on the selection band, and 
did not characterize uniquely a source. Recent all-sky X-ray 
surveys (Laurent-Muehleisen et al. 1993; Brinkmann et al. 
1995; Perlman et al. 1995a), in fact, include many BL Lacs 
previously selected at radio frequencies, which could now 
then be classified both as XBLs or RBLs. 

It has been known for some time (e.g., Ledden & O'Dell 
1985) that the broad-band spectra of the two classes are also 
different (with the likely exception of the radio to millime- 
tre band: Gear 1993). Recently, Giommi, Ansari & Micol 
(1995a) have studied the multifrequency spectrum of a large 
sample of BL Lacs, arguing that there may be only one pop- 
ulation of objects, characterized by a wide range of energy 
cutoffs (or alternatively energy peaks in a ff^ vs. v plot): at 
infrared-optical energies for most RBLs, at ultraviolet/X-ray 
or higher energies for most XBLs, with the latter constitut- 
ing only ~ 10 per cent of all BL Lacs. The X-ray-to-radio 
flux ratios (/x//r) for the two classes are also different, as 
a result of the different break frequencies. Note that within 
this hypothesis there is a continuum of cutoff energies, with 
different selection bands picking up objects with different 
spectra. 

We (Giommi & Padovani 1994; Padovani & Giommi 
1995a) have then proposed a division based on /x//r (as the 
break frequency at present can be derived directly from the 
spectra only for a small number of objects: see Sect. 3), clas- 
sifying as HBLs (or high-energy cutoff BL Lacs) and LBLs 
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(or low-energy cutoff BL Lacs) sources respectively above or 
below the value of /x//r ~ 10^^^ (where X-ray fluxes cover 
the 0.3 - 3.5 keV range and are in units of erg cm~^ and 
radio fluxes refer to 5 GHz and are expressed in janskys). 
We have also shown that the hypothesis that HBLs (that 
is, so far, most XBLs) are the minor component of the BL 
Lac population with high-energy cutoffs can explain most of 
their properties at least as well as the alternative (and more 
diffuse) explanation for the existence of two classes of BL 
Lacs, which posits that XBLs arc intrinsically more numer- 
ous as they arc viewed significantly off their beaming axis 
(e.g., Stocke ct al. 1985; Maraschi et al. 1986; Padovani & 
Urry 1990). 

A study of the puzzling distribution of BL Lacs on the 

Qro,aox plane (Padovani & Giommi 1995a) has led us to 
suggest that the observed bimodality can also be explained 
by an energy break ranging from high (X-ray) to low (far- 
infraxed) energies for HBLs and LBLs respectively. This as- 
sumes a different origin of the X-ray emission for HBLs 
and LBLs: namely, an extension of the synchrotron emission 
likely responsible for the lower energy continuum in the for- 
mer and inverse Compton (SSC: synchrotron self-Compton) 
emission in the latter. This would be consistent with the 
peak of the emitted power being at lower energies for LBLs, 
suggestive of an extra component taking over at higher fre- 
quencies (see Bregman 1990 and references therein). A cru- 
cial test of our hypothesis (and in general of amy other model 
wishing to explain the existence of two BL Lac classes) is 
then the analysis of the spectral shape of the X-ray emission 
of a large sample of BL Lacs including both HBLs and LBLs. 
In the X-ray band, in fact, Compton emission is expected to 
be flatter than the synchrotron component (e.g., Ghisellini 
& Maraschi 1989). 

The study of X-ray spectral indices of BL Lacs has given 
various and sometimes contradicting results, even consider- 
ing only the most recent papers dealing with possible differ- 
ences between the two classes (references to earlier papers 
can be found in Sambruna et al. 1994). Sambruna et al. 
(1994) have analyzed the EXOSAT spectra of 21 BL Lacs 
(only five of which were RBLs) and found similar spectral 
indices in the LE (0.04 - 2.0 keV) and ME (1 - 20 keV) 
bands for the two classes, with indications of a steepening 
at higher energies for XBLs. Ciliegi, Bassani & Caroli (1995) 
have studied the X-ray spectra of 42 BL Lacs (mostly RBLs) 
observed with various {pve-ROSAT) X-ray telescopes in dif- 
ferent energy ranges and found that XBLs and RBLs have 
similar (within the rather large errors) spectral indices in the 
soft (0.2 — 4 keV) and possibly also in the hard (2 — 10 keV) 
band. However, there were indications that XBLs and RBLs 
have different global X-ray spectra, with the former and the 
latter steepening and flattening, respectively, at higher en- 
ergies. Urry et al. (1996; see also Sambruna 1994) have ex- 
amined more recent _R05Ar spectral data for 32 out of the 
34 radio-selected BL Lacs in the 1 Jy sample (Stickel et al. 
1991). They flnd that 1 Jy RBLs have spectra similar to 
those of the Einstein Observatory Extended Medium Sen- 
sitivity Survey (EMSS) XBLs (Perlman et al. 1995b), with 
an indication for the subsample with Einstein IPC data of 
a flattening of the spectra at higher (> 1 keV) energies. Co- 
mastri, Molendi & Ghiselhni (1995) linaly zed the ROSAT 
spectra of twelve 1 Jy BL Lacs (i.e. a subsample of the 
objects studied by Urry et al. 1996) dividing the objects 



according to the X-ray-to-radio flux ratio, as proposed by 
Padovani & Giommi (1995a), with the result that the aver- 
age spectrum of the (3) HBLs was steeper than that of the 
(9) LBLs. Urry et al. (1996) reached similar conclusions. 

It is clear that at least some of the different results can 
be explained by poor statistics and the use of the RBL/XBL 
classification, which as we have discussed above is not phys- 
ical. However, it is also probably true that the differences 
between the X-ray spectra of HBLs and LBLs are not that 
strong and require better statistics to be singled out. 

The purpose of this paper is to analyze the X-ray spec- 
tra of all BL Lacs observed (as pointed or serendipitous 
sources) by ROSAT to check for any differences between 
HBLs and LBLs, coupling for the first time a large statis- 
tics with what we think is a meaningful distinction between 
the two classes. Our data have been taken from the WGA 
catalogue (White, Giommi & Angelini 1994), a large list of 
X-ray sources generated from all the ROSATVSPC pointed 
observations, with which is now possible to study the X-ray 
properties of large numbers of objects in an homogeneous 
and relatively simple way. The selection of the objects was 
done by cross-correlating the first revision of the WGA catar 
logue (restricting ourselves to sources with quality flag > 5, 
which excludes problematic detections) with our recent cat- 
alogue of BL Lacs (Padovani & Giommi 1995b). This pro- 
duced 163 observations of 85 distinct BL Lacs (58 HBLs 
and 27 LBLs: see Sect. 2), which correspond to ~ 46 per 
cent of confirmed BL Lacs presently known. This represents 
the largest number of BL Lacs for which homogeneous X- 
ray spectral information is available and the largest BL Lac 
sample ever studied at X-ray frequencies (cf. Ciliegi et al. 
1993 who collected information on the spectra of 42 BL Lacs 
observed with various X-ray telescopes). 

The structure of the paper is as follows: Section 2 de- 
scribes the observational data and the data analysis. Sec- 
tion 3 studies the X-ray spectral properties of BL Lacs 
while Section 4 discusses our results and presents our con- 
clusions. Throughout this paper spectral indices are written 



2 DATA ANALYSIS 

Spectral indices for the WGA BL Lacs were obtained from 
the count rates given in the catalogue in the 0.1 — 0.4 keV 
range (soft band: S), 0.4 — 0.86 keV range (mid band: M) 
and 0.87 — 2.0 keV range (hard band: H) combined into 
hardness (HR) and softness {SR) ratios (see Giommi et 
al. 1995b). Briefiy, the count rates were combined to con- 
struct one SR = S/M and two HRs, HRi = H/M and 
HR2 = H/{M + S). Initially we converted the hardness ra- 
tios into energy spectral indices both assuming Galactic A'h 
derived from 21 cm measurements (Stark et al. 1992; Shafer 
et al., private communication) and with A^h derived from the 
softness ratio through an iterative procedure. The spectral 
indices derived from the two approaches were found to be 
similar and well correlated in the 0.4 — 2.0 keV range, while 
considerable scatter was present when the total ROSAT 
band was considered, with a sizeable number of objects hav- 
ing ax,NH significantly different from Ox.Gaiactic nh, sugges- 
tive of a change of the spectrum at lower energies. In this pa- 
per we will adopt as X-ray spectral index the value obtained 
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with Nyi fixed to the Galactic value, which is better deter- 
mined than the A''h derived from the softness ratios, and in 
the 0.4 — 2.0 keV range, i.e. the mid to hard band (spectral 
indices in the whole ROSAT range were used for the com- 
parison between the values obtained with our method and 
those derived by other authors from detailed spectral fits: 
see below). 

The choice of this energy range stems from two facts: 
first, our method assumes a single power-law; given the ev- 
idence of a change of spectral shape at lower energies in 
some sources, as described above (see also Comastri et al. 
1995 and Urry et al. 1996), this assumption might be incor- 
rect in the whole ROSAT range while it is probably more 
acceptable in the narrower 0.4 — 2.0 keV range (as shown by 
the similarity of the spectral indices obtained with Galac- 
tic A^H and A^h derived from the softness ratio) ; second, the 
method used in the WGA catalogue to estimate the source 
intensity uses the counts detected in a box whose size opti- 
mises the signal to noise (S/N) ratio. This size is calculated 
assuming an average point spread function (PSF) that is too 
sharp for very soft photons. For weak sources near the field 
center the "optimum" box size is small and consequently 
includes a small, energy-dependent, fraction of the source 
photons causing an underestimation of the counts in the soft 
band. The effects of the dependence of energy response with 
off-ajcis angle have been taken into account using different 
conversion matrices in five different off-axis ranges: — 20, 
20 - 30, 30 - 40, 40 - 50, and 50 - 60 arcminutes. 

Errors on the spectral indices (Icr) were derived from 
the uncertainties on the hardness ratios. We included in our 
analysis only observations with values of hardness ratios rel- 
ative to their uncertainties > 3.5, which corresponds roughly 
to 1(7 errors on fix <; 0.5 (Giommi et al. 1995b). All objects 
in this study arc "bona fide" BL Lacs: we excluded eight 
sources classified as uncertain BL Lacs or BL Lac candi- 
dates (Padovani & Giommi 1995b). We are then left with 
163 observations of 85 BL Lacs. 

Figure 1 compares the X-ray energy spectral indices 
derived from our method in the whole ROSAT band (i.e. 
obtained from HR2) with those obtained by other authors 
(Perlman et al. 1995b; Maraschi et al. 1995; Urry et al. 1996) 
through a detailed spectral fitting procedure, using the same 
ROSAT data for 44 BL Lacs. The comparison is between 
spectral indices derived under the same assumption regard- 
ing the absorbing column density, i.e. either assuming Galac- 
tic A'h (Stark et al. 1992; Shafer et al., private communica- 
tion) or with A'h derived from the softness ratio. Figure 1 
shows that most objects cluster around the line Qx.hr = Qx.flt 
with a linear relationship between the spectral indices de- 
rived with the two methods. Overall the agreement is then 
very good, especially considering the underestimate of soft 
photons in the WGA data. The mean values for the two 
derivations are Ox.hr = 1.18 ± 0.08 and Ox.flt = 1.24 ± 0.08 
(here and in the following we give the error of the mean), 
which clearly indicates the statistical robustness of our ap- 
proach. 

A few sources in the PSPC field of famous objects have 
many X-ray observations: for example, IE 1415+2557 (in 
the field of NGC 5548) has been observed in 29 different 
occasions! Inclusion of all these data in statistical studies 
(means, correlations, Kolmogorov-Smirnov [KS] tests etc.) 
would clearly bias them: the resulting mean, for example. 
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Figure 1. ROSAT band energy spectral indices of BL Lacs de- 
rived from our approach based on the hardness ratios, Qhr; com- 
pared to the values obtained by other authors from a detailed 

spectral fitting procedure of the same ROSAT data, af\t. The 
solid line indicates the locus of points with ajj^ = Qflt . Error bars 
represent la errors for our values and 90 per cent errors for those 
obtained by other authors. 

would represent more the typical X-ray spectrum of IE 
1415-1-2557 than that of BL Lacs as a class. Therefore, al- 
though we plot in the figures all the 163 X-ray spectral in- 
dices, in all statistical tests we have kept only one observs/- 
tion (and therefore one spectral index) per object, selected 
on a one-to-one basis taking the best combination of offset 
from the field center and S/N ratio. When more than two 
observations with comparable offset and S/N ratios were 
available, we selected as most representative the one with 
ax nearest to the mean value for that source. 



3 THE X-RAY SPECTRA 

3.1 Ox versus /x//r 

The main purpose of this paper is to study the X-ray spec- 
tral indices of HBLs and LBLs. To this aim, it is important 
to define the value of the /x//r ratio which divides the two 
classes. (We prefer the use of broad-band X-ray fluxes in 
the X-ray-to-radio flux ratios instead of, for example, 1 keV 
fluxes, because the latter are more dependent on the precise 
value of Qfx and therefore require spectral fitting, which is at 
present not possible for all BL Lacs.) In Padovani & Giommi 
(1995a), based on the energy distributions given by Giommi 
et al. (1995a), we suggested a ratio /x//r ~ 10~^^. It should 
be noted, though, that very few BL Lacs in the compilation 
by Giommi et al. (1995a) had lO"^^'''' < /x//r < 10""' ■^ be- 
cause this intermediate zone is not favoured either by radio 
or X-ray selection and is therefore poorly populated both in 
X-ray and radio surveys. The dividing /x//r value was then 
not well determined. 
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Figure 2. The X-ray spectral index of BL Lacs versus their 
X-ray-to-radio flux ratio (with X-ray fluxes in the 0.3 - 3.5 keV 
range in units of erg cm~^ and radio fluxes at 5 GHz ex- 
pressed in janskys). Error bars represent Icr errors. 

In Fig. 2 we plot Ox versus /x//r for all the objects in our 
sample, with radio data from Padovani & Giommi (1995b 
and references therein) and X-ray fluxes in the 0.3 3.5 keV 
band derived from the Galactic A'h and values. (Ide- 
ally the radio and X-ray observations should be simultane- 
ous; although variability in the radio is not very strong, this 
undoubtedly introduces some "noise" in our plots.) Within 
our hypothesis on the difference between HBLs and LBLs 
(Padovani & Giommi 1995a), described in the Introduction, 
the X-ray emission of HBLs should be dominated by syn- 
chrotron emission. HBLs near the transition region, in which 
the break is just below the X-ray band and where we see the 
steep tail of the synchrotron emission but not yet the inverse 
Compton component, should then be characterized by the 
steepest X-ray spectra (see also Sect. 3.3). Figure 2 shows 
this to happen for lO""'^ < /x//r < 10"", where aU ob- 
jects have Ox > 1.5 and reach values up to ~ 3. A better 
dividing line between the two classes is then /x//r ~ 10~^^'^ , 
which is the value at which there is a discontinuity in Qx- 
We stress that our results (and those obtained in Padovani 
& Giommi 1995a) are unaffected by the precise value of this 
parameter. (One object has ax ~ 3 but /x//r ~ 10~^^. This 
is IWGA J1202. 1-1-444, which appears three times in the 
WGA catalogue and only in its faintest state crosses the 
HBL/LBL boundary to become an LBL with an unusually 
steep X-ray spectrum, albeit with large errors, following the 
Qx — /x anti-correlation normally observed in HBLs [e.g., 
Giommi et al. 1990; Sambruna et al. 1994]. We note how- 
ever that a decrease in radio flux by a factor of 2 — 3 as 
compared to the value we used, which we camiot rule out, 
would have kept this object within the HBL limits.) 

The mean value of the X-ray spectral index for the 85 
BL Lacs under study is (a^) = 1.37 ± 0.05. Figure 2, how- 
ever, already shows that the X-ray spectral indices of the 



58 HBLs (BL Lacs with /x//r > 10"" '^) are steeper than 
those of the 27 LBLs (BL Lacs with /x//r < lO"""""'). In fact, 
(ox.hbl) = 1.52±0.06 while (qx.lbl) = 1.06±0.09, different 
at the 99.99 per cent level according to a Student's t-test. 
The two ax distributions, shown in Fig. 3, axe also differ- 
ent at the 99.9 per cent level according to a KS test. (Had 
the dividing line been drawn at the value /x//r = 10"'^'^, 
against the evidence presented in Fig. 2, one would have 
found (ax.HBL) = 1.43 ± 0.05 and (ax.LBL) = 1.27 ± 0.11, 
with the two distributions different at the 95.4 per cent level 
according to a KS test.) 

Another notable feature of Fig. 2 is the anti-correlation 
between ax and /x//r for HBLs, significant at the 99.99 per 
cent level (Spearman rank-order test), independently of the 
precise value of the dividing /x//r ratio, and not present in 
LBLs. A linear fit to the HBL data gives Qx = —(0.41 ± 

0. 07) Iog(/x//r) - (2.65 ± 0.66). In HBLs, in other words, 
the smaller the X-ray to radio fiux ratio, (i.e., the steeper 
the effective radio- X-ray spectral index), the steeper is the 
X-ray spectral index, which points to a connection between 
the radio and X-ray bands and/or to the effects of a break 
in the energy distribution (see Sect. 3.3). 

3.2 Ox versus Qox 

A similar and related connection for HBLs exists between 
the optical and X-ray bands. Figure 4 plots the X-ray spec- 
tral slope versus aox, the effective optical- X-ray spectral in- 
dex evaluated between the rest-frame frequencies of 5000 
A and 1 keV. Optical fiuxes were derived from the V mag- 
nitudes collected by Padovani & Giommi (1995b), correct- 
ing for Galactic absorption following Giommi et al. (1995a), 
while 1 keV fiuxes have been obtained from the ROSAT 
counts and the derived spectral indices and have also been 
corrected for Galactic absorption. The fc-correction has been 
derived assuming an optical index ao = 1.05 for LBLs and 
ao = 0.65 for HBLs (Falomo, Scarpa & Bersanelli 1994). For 
fc-correction purposes objects with no redshift have been as- 
signed the mean redshift of the corresponding sub-sample, 

1. e. amean.LBL = 0.60 and 2mean,HBL = 0.25. Notc that Opti- 
cal fiuxes are not simultaneous with X-ray data (although 
for most of the 1 Jy BL Lacs in our sample they should be 
representative of their "typical" optical state: see Padovani 
& Giommi 1995b), which will certainly introduce a scatter: a 
variation AV in magnitude, in fact, translates into a change 
Aaox = 0.15AV. 

From the data presented in Fig. 4 we infer that: 1. for 
HBLs the X-ray spectral index is correlated with Qqx at the 
99.99 per cent level (Spearman rank-order test) with a best 
fit ax = (1.38 ± 0.20)aox - (0.01 ± 0.22); 2. no correlation 
is present for LBLs; 3. the two classes separate extremely 
well on this plane, that is all HBLs apart from one have 
Ox > Oox (within the errors), while the majority of LBLs 
have Qx < Qox- 

This separation is better appreciated in Fig. 5, which 
plots Ox — ctox versus fx/ fi- This plane can be divided into 
four quadrants, out of which one is basically unpopulated 
and another one is only marginally populated. The remain- 
ing ones are the HBL-convex-spectrum quadrant {fx/ fv > 
10~^^'^ and Ox > aox) and the LBL-concave-spectrum 
quadrant (/x//r < lO""'^ and Qx < Qox). The only HBL 
with possibly concave spectrum is MS0122. 1-1-0903, which 
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Figure 3. The X-ray spectral index distribution of all BL Lacs 
m our sample (dotted line), and of HBLs {/x//r > lO-"-^: solid 
line) and LBLs (/x//r < dashed line) separately. Only 

one X-ray spectrum per object is used, as described in Sect. 2. 
The distributions of HBLs and LBLs are different at the 99.9 per 
cent level according to a KS test. 
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Figure 4. The X-ray spectral index versus ctox for BL Lacs. Filled 
points indicate HBLs (/x//r > 10~^^'^), while open squares indi- 
cate LBLs (/x//r < 10~^^'®). The solid line represents the locus 
of points having Qx = Qox- Error bars represent Itr errors. 

is however consistent with having ax > ccox at the ~ 20- 
level, especially considering that errors on ax — aox do not 
include errors on aox- As regards the twelve LBL spectra 

with possibly convex optical- X-ray continuum, only four of 
them have Qx > Oox at the > 2cj level. These correspond to 



four different objects: one is IWGA J1202. 1-1-444, discussed 
above, while the other three have /x//r values within a fac- 
tor of two of the HBL/LBL dividing line, which, as far as 
we know, could have well been crossed had there been a de- 
crease in radio flux by a factor of 2 (compared to the values 
we used) at the time of the ROSAT observation. 

The difference between HBLs and LBLs in (qx — Qox) is 
striking, this parameter being equal to 0.40 ± 0.04 for HBLs 
and -0.33d=0.09 for LBLs (0.17±0.06 for the whole sample). 
This is simply an amplification of the different X-ray and 
effective optical- X-ray spectral indices. While, as we have 
seen, HBLs have steeper than LBLs, the opposite is true 
as regards aox (see Fig. 4), for which (aox.HBL) = 1.11±0.03 
and (aox.LBL) = 1-39 ± 0.03. 

Sambruna (1994; sec also Maraschi et al. 1995) has 
found a significant difference between the Qox — Qx distri- 
butions of 1 Jy RBLs and EMSS XBLs, with mean values, 
however, apparently less different than ours. This difference 
can be explained by her smaller statistics and inclusion of 
three 1 Jy HBLs with the RBL sample. 

Perlman et al. (1995b) have compared the ROSAT X- 
ray spectral indices of 22 EMSS BL Lacs (all HBLs) to the 
slope of the tangent at 1 keV to a parabola fitted to the 
radio, optical, and X-ray fluxes. They find that the X-ray 
spectra have similar or slightly steeper slopes than the tan- 
gent slope in most cases, a somewhat different result from 
ours. This might be due to our larger sample and also to their 
use of the tangent slope instead of the effective opticaJ-X- 
ray index. Brunner et al. (1994) have studied the ROSAT 
spectra of five 1 Jy BL Lacs and derived a mean value for 
Ox — Oox consistent with zero, unlike our result for LBLs. 
Despite their small statistics, exclusion of the HBL object 
in their sample (S5 0716 + 714), which has the largest value 
of Ox — Oox, makes their mean value consistent with ours. 

3.3 ax versus i^ireak 

The ant i- correlation between Qx and /x//r and the correla- 
tion between Ox and Oox for HBLs are both related to what 
we think is the most interesting correlation, that between 

ax and inbreak- 

In theory, one should infer the break frequency i^ieak 
from the observed energy distribution of the objects. In prac- 
tice, however, the available spectral coverage is far from com- 
plete and the measurements are generally non-simultaneous. 
In particular, far-ultraviolet and far-infrared data are miss- 
ing for the large majority of the sources, and these are the 
bands where the peak of the emission is most likely located 
for HBLs and LBLs respectively (Giommi et al. 1995a). The 
brealc frequency has then been derived as follows: Padovani 
& Giommi (1995a), using a triple power-law approximation 
to the spectra, have shown that the distribution of BL Lacs 
on the Qro, Qox plane (where Qro is the effective radio-optical 
spectral index evaluated between the rest-frame frequencies 
of 5 GHz and 5000 A) can be easily explained by an energy 
break migrating from high (X-ray) to low (far-infrared) en- 
ergies (see their Fig. 12). Here we use their results to infer 
inbreak from the Qro and Qox values. The former were derived 
from the radio and optical data collected by Padovani & 
Giommi (1995a), using optical indices and mean redshifts 
(when necessary) as described above. As regards the radio 
spectral index, necessary for the fc-correction of Qro, for the 
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Figure 5. ax minus Qox versus /x//r for BL Lacs. The horizontal 
dashed line represents the locus of points with Ox = aox, while 
the vertical dashed line represents the dividing line between LBLs 
and HBLs at /x//r = 10^^^'®. Error bars represent la errors and 
include only the uncertainties on ax. 

BL Lacs lacking this parameter it was assumed equal to the 
mean value for the two sub-classes, that is Or = 0.2 for HBLs 
and Or = —0.2 for LBLs. 

Operationally, we projected the points on the 
plane on the predicted relationship between the two param- 
eters derived by Padovani & Giommi (1995a), namely on 
that of HBLs for /x//r > lO-" '^ and on that of LBLs for 
/x//r < 10-"-^ (see Fig. 12 of Padovani & Giommi 1995a). 
We then read off the corresponding values of the break fre- 
quency. Although this method gives only an approximate 
estimate of Vbieak, the values derived this way agree well 
with those inferred from an inspection of the plots, for ex- 
ample, in Giommi et al. (1995a) for the (few) objects with 
enough data to permit a direct determination of the break. 
Moreover, our values are typically consistent within ~ 0.2 
dex with the peak frequencies obtained by Comastri et al. 
(1995) from polynomial fits to the broad band spectra of 
twelve 1 Jy BL Lacs. 

Figure 6 clearly shows the strong anti-correlation (sig- 
nificant at the 99.99 per cent level: Spearman rank-order 
test) between X-ray spectral index and break frequency for 
HBLs. A linear fit gives Qx = —(0.63 ± 0.09) log inbreak -|- 
(11.87 ± 1.50): i.e., the highest inbreak, the flattest Ox. As 
regards LBLs, here the situation is not clear: a fit to the 
"compressed" data as done for all the correlations studied in 
this paper (i.e., using only one spectrum per object, selected 
as described in Sect. 2) gives no correlation between Qx and 
inbreak for thc 27 LBLs (Pcorr ~ 61 per cent). However, 
Fig. 6 shows an apparent positive correlation between the 
two parameters (qx oc 0.411ogi'break) when all (42) spectra 
are included, significant at the 99.6 per cent level (Spear- 
man rank-order test; this reduces to 98.2 per cent, with 
Qx oc 0.26 log I'break, if One excludes the two outliers with 
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Figure 6. The X-ray spectral index versus the logarithm of the 

break frequency for BL Lacs. Objects with inbreak > 10^'"' Hz have 
/x//r > 10~^^'^ and are, by definition, HBLs, while the remaining 
objects arc LBLs. Error bars represent Icr errors. See text for the 
derivation of t-break- 

thc largest and smallest spectral index). This might suggest 
that LBLs display a (weak) positive correlation, possibly di- 
luted by the relatively small statistics and larger error bars 
(the LBL sample is less than half the size the HBL one and 
on average the LBLs studied in this work are weaker X-ray 
sources than HBLs). 

3.4 Interpretation 

The negative correlation between Qx and i^broak observed for 
HBLs is expected if the emission of these objects is due to 
the synchrotron process all the way up to the X-ray band, as 
shown in Fig. 7: when the break is near the X-ray band (~ 
10^^ Hz: dotted line) thc X -ray spectrum will be relatively 
flat, but as soon as thc peak shifts to lower energies (~ 10^® 
Hz: dashed line) the steep tail of the synchrotron emission 
moves in. When fbreak is at much lower energies (e.g., ~ 
lO'^* Hz: solid line) so that the synchrotron component is no 
more dominant, Compton omission starts to be important: 
at this point one has crossed the HBL/LBL dividing line. For 
even smaller break frequencies, two effects compete: first, the 
smaller lAjieak, the smaller is the product S7max,ej where B is 
the magnetic field and 7max,e is the maximum Lorentz factor 
of the electrons (e.g., Rybicki & Lightman 1979); this should 
not affect the lower energies, thereby leaving unchanged the 
shape of the Compton component; second, a lower ;4,reak 
means a weaker steep X-ray synchrotron component, which 
should result in a flatter spectrum in objects with still a non- 
negligible synchrotron omission. The sum of the two effects 
is not easy to quantify but the fact that in LBLs the X-ray 
spectral index is not strongly dependent on z/break suggests 
that the Compton component might dominate completely 
the X-ray band in at least some objects. 
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This ties in nicely with the study of the Qox index, which 
shows that the X-ray band connects smoothly to the opti- 
cal one in HBLs, with a steepening of the broad-band spec- 
trum, supporting a common emission mechanism for the two 
bands, while in the great majority of LBLs Ox is flatter than 
Oox, indicative of a different component entering the X-ray 
range. Also, the anti-correlation between Qx and /x//r in 
HBLs is simply related to that between and t'break as 
objects with larger inbreak have also larger values of f^/ fr- 

In objects with very low break frequencies, for which the 
X-ray emission should be pure inverse Compton, the X-ray 
spectral indices should be similar to those in the radio-far- 
infrared region, obviously up to f < inbreak- Unfortunately, 
very few BL Lacs have been detected by IRAS so the statis- 
tics is quite poor: we find that the 5 GHz - 60/i spectral index 
for about ten 1 Jy BL Lacs is Ori ~ 0.3 (sec also Giommi 
et al. 1995a), which roughly corresponds to the lower bound 
of Qx for LBLs. Clearly a detailed study of the far infrared 
spectra of LBLs would be extremely valuable in this respect: 
this project will be carried out with ISO. 

We predict that the observed anti-correlation between 
Qx and I'brcak obscrvcd for HBLs cannot continue indefi- 
nitely. Since the synchrotron spectrum steepens with fre- 
quency (e.g., Ghisellini & Maraschi 1989), if Vbreak moves to 
larger and larger values, the X-ray spectral index at min- 
imum will be as flat as the effective radio-optical spectral 
index Oro- In that case, then, one will have ax ~ Qox ~ 
Oro ~ 0.3 — 0.5, the typical values for HBLs. Extrapolating 
the observed anti-correlation, one should find Ox ~ 0.4 ~ 
constant for t'brcak J> 10 kcV. As this requires very high 
electron energies, it is not surprising that no BL Lac in our 
sample has such a high I'break- If such objects exist they are 
probably quite rare but should be detected by future hard- 
X-ray missions. 

We note that the paucity of objects with break frequen- 
cies in the optical/ultraviolet band lO"'^ - lO^'^'^ Hz) is 
simply a selection effect, as these objects have intermediate 
values of /x//r which make them difficult to detect both in 
radio and X-ray surveys. Optical surveys, which in princi- 
ple do not discriminate against this sort of objects, are at 
present very inefficient (e.g.. Green, Schmidt & Liebert 1986; 
Jannuzi, Green & French 1993) and have so far led to the 
selection of only about ten BL Lacs (Padovani & Giommi 
1995b). 

4 DISCUSSION AND CONCLUSIONS 

The data presented in this paper allow us to draw a con- 
sistent picture of HBLs and LBLs. Namely, the correlations 
of the X-ray spectral indices with Oox, /x//r, and break fre- 
quency all point to a common, likely synchrotron, origin 
for the radio/optical/X-ray emission in HBLs, with an over- 
all convex spectrum. On the other hand, as hinted by the 
lower cutoff energies, in LBLs a second, harder, most likely 
inverse Compton, component appears in the X-ray band. 
The ROSAT spectral indices of the two classes are differ- 
ent, with HBLs displaying on average steeper spectra which 
show a strong anti-correlation with break frequency, while 
LBL spectra might show a positive dependence on this par 
rameter. These are yet other differences between the two 
BL Lac classes, which add to the previously known ones, 
and might be helpful in selecting BL Lac candidates from 
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Figure 7. Examples of multifrequency spectra of BL Lacs with 
different peak (or break) frequencies: ~ lO^'^ Hz (~ 0.4 keV: 
dotted line), ~ IQl^ Hz (dashed line), ~ lO" Hz (solid line). 
The two parallel (dot-daslicd) lines represent the ROSAT ha,nd. 
Note how the X-ray spectrum steepens with decreasing break 
frequency until Compton emission becomes the dominant process. 
The normalizations of the different curves have been chosen to 
best illustrate the spectral differences in the X-ray band and do 
not reflect intrinsic properties. 

the WGA catalogue, for example. 

Based on the additional information presented here, we 
can better pinpoint the borderline between HBLs and LBLs, 
and therefore propose a refined dividing line between the two 
classes at /x//r ~ 10~^^'^, which corresponds to arx ~ 0.8 
between 1 keV and 5 GHz and a break frequency ~ 10^^ Hz 
(see Fig. 6). Our results agree with the hypothesis we put 
forward that the cutoff energy is the fundamental parameter 
giving rise to two different classes of BL Lacs (Giommi & 
Padovani 1994; Padovani & Giommi 1995). Oddly enough, 
our results are also consistent with the other explanation 
for the existence of two BL Lac classes based on orientation. 
Ghisellini & Maraschi (1989), in fact, have shown that in 
objects at larger angles with respect to the line of sight (i.e., 
XBLs/HBLs in their view) synchrotron emission should be 
the main radiation mechanism in the X-ray band, while in- 
verse Compton should be more important for more aligned 
sources (i.e., RBLs/LBLs), in agreement with our findings. 
The main problem of that model, however, apart from other 
ones discussed by Padovani & Giommi (1995a), is that a 
change in viewing angle does not seem to be able to ex- 
plain the large difference in cutoff energies between the two 
classes, as exemplified by Fig. 4 of Ghisellini & Maraschi 
(1989) (Sambruna 1994; Sambruna et al. 1995, in prepara- 
tion). 

Our results explain the peculiar spectral variability of 
XBLs (i.e., HBLs) and RBLs (i.e., LBLs). Various studies 
have shown that in XBLs higher fiux states are associated 
with flatter spectra (e.g., Giommi et al. 1990; Sambruna et 
al. 1994). Little is known about spectral variability in RBLs 



8 P. Padovani and P. Giommi 



but Urry et al. (1996) have shown that the handful of 1 
Jy BL Lacs with suitable data seem to show the opposite 
trend, that is steeper spectra in brighter states (the excep- 
tion, PKS 2005 - 489, is an HBL). If X-ray variability in BL 
Lacs is accompanied by a change in brealc frequency (e.g., 
George, Warwick & Bromage 1988: MKN 421; Sambruna 
1994: S5 0716+714, PKS 2005-489), this behaviour is eas- 
ily explained: when an object brightens and the break moves 
to higher energies, in HBLs the spectrum hardens, as shown 
in Fig. 8 for a single source and in Figs. 2 and 6 for the HBL 
population. The opposite effect for RBLs can be explained 
in those cases where synchrotron emission in the X-ray band 
is not negligible in the high state: when the break moves to 
higher energies, the steep synchrotron component moves in 
the X-ray band therefore steepening the spectrum (Fig. 8). 
This would be in agreement with the (weak) correlation be- 
tween Ox and Inbreak for LBLs, discussed in Sect. 3.3 (see Fig. 
6) . Figure 8 also shows that a spectral variability similar to 
that observed in the X-ray band for HBLs, that is flatter 
spectra at higher states, should be mirrored by LBLs in the 
infrared band, as in fact observed (e.g.. Gear ct al. 1986). 

Our results seem in contradiction with those of Urry et 
al. (1996), who found that the ax distribution for 1 Jy RBLs 
was consistent with that of the EMSS XBLs. As discussed in 
the Introduction, however, most previous work on the study 
of spectral indices of BL Lacs was limited by smaller statis- 
tics and the use of the RBL/XBL classification, which can 
be misleading. In particular, the result of Urry et al. (1996) 
can be explained as follows: 1. three of the 1 Jy RBLs are 
HBLs (S5 0716 + 714, MKN 501, PKS 2005 - 489), whose 
inclusion in the RBL sample certainly shifted the Ox dis- 
tribution to steeper values (the mean value of for these 
three objects is significantly different, i.e., steeper, than that 
of the rest of the sample: see also Comastri ct al. 1995 and 
Urry et al. 1996); 2. the EMSS BL Lacs (all HBLs) have val- 
ues of /x//r ^ 10-^°-^ (Padovani & Giommi 1995a), partly 
because, being a serendipitous survey, the EMSS excludes 
by definition well-know BL Lacs with considerable radio 
emission which have lower /x//r values but more impor- 
tantly because objects with relatively low values of /x//r 
have a small probability of being detected and large-area 
X-ray surveys are needed. (This interpretation is supported 
by the fact that the Slew survey HBLs reach lower /x/Jr 
values: see Perlman et al. 1995a). Given the anti-correlation 
between and /x//r (Fig. 2), this implies that EMSS HBLs 
are bound to have ax flatter than average, thereby reducing 
any difference between the two classes. In fact, if we com- 
pare our LBL sample with the EMSS HBLs, wo also find no 
significant difference between the spectral indices of the two 
samples. 

One could argue that our sample is not complete in a 
statistical sense. However, it is the largest sample for which 
homogeneous X-ray data are available and includes 50 per 
cent of known BL Lacs. We needed a large enough sample 
to maximize the coverage of the parameter space to be able 
to test for any differences in the X-ray spectra of the two 
BL Lac classes. Previous results, even those based on com- 
plete (and therefore not very large) samples, were in fact not 
clear-cut. We are not aware of selection effects which could 
possibly alter our results significantly. Furthermore, consid- 
ering only the complete sub-sample of 22 EMSS BL Lacs 
(all HBLs: Perlman et al. 1995b) we still reproduce one of 
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Figure 8. The effect of a change in break frequency on the flux 
and spectral variability of BL Lacs: in HBLs (dashed lines), when 
the object brightens and the break moves to higher energies, the 
X-ray spectrum hardens. In LBLs (solid lines), the opposite ef- 
fect might apply (i.e., spectral steepening at higher fluxes) wiien 
synchrotron emission in the X-ray band is not negligible in the 
high state. Incidentally, note also the much larger flux variability, 
at fixed radio variability, expected for LBLs in the optical band 
(parallel dotted lines), as compared to that of HBLs, as effectively 
observed (e.g., Jannuzi et al. 1994). As in Fig. 7, the normaliza- 
tions of the different curves have been chosen only for illustrative 
purposes and do not reflect intrinsic properties. 

our most important results, i.e. we find an anti-correlation 
between Qx and z^broak, significant at the 99.3 per cent level 
(Spearman rank-order test), consistent with the one we ob- 
tained for the whole HBL sample. 

Comastri et al. (1995) have found a positive correlation 
between Qx and inbreak in their sample of twelve 1 Jy BL Lacs 
(which includes nine LBLs and three HBLs). This is consis- 
tent with our results. If we consider all objects irrespective 
of the HBL/LBL classification, in fact, we also detect a pos- 
itive correlation between Ox and inbreak, because on average 
HBLs have steeper spectra. Our much better statistics, how- 
ever, allows us to examine this relation in more detail and 
to see that the dependence of Ox on break frequency is dif- 
ferent for the two BL Lac classes, in particular HBLs exhibit 
a negative correlation between the two parameters. 

The main conclusions of this paper, which studies the 
ROSATXrvecy spectra of about half the known BL Lacs, can 
be thus summarized: 

BL Lacs are characterized by energy power-law spectral 

indices ranging between ~ and 3, with an average value 
~ 1.4. The X-ray spectral properties of BL Lacs with high 
(> lO^-""' Hz) frequency cutoffs (HBLs) are significantly differ- 
ent from those with lower frequency cutoffs (LBLs). HBLs 
have on average steeper X-ray spectral indices (ox ~ 1.5) 
which correlates with the effective radio-X-ray and optical- 
X-ray spectral indices. The overall spectrum is convex from 
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radio to X-ray energies, and gets steeper with frequency. A 
strong anti-correlation is also present between X-ray spec- 
tral index and break frequency, with sources having cutoffs 
at larger energies displaying harder X-ray spectra. In LBLs, 
on the other hand, which have ax ~ 1.1, the X-ray spec- 
trum shows no correlation with the effective spectral indices 
and the optical-X-ray continuum has a concave shape. These 
differences might also explain the peculiar X-ray spectral 
variability of HBLs and LBLs. 

These findings strongly support the hypothesis that the 
main radiation mechanism in HBLs is synchrotron emission 
from the radio to the X-ray band, while in LBLs an inverse 
Compton component dominates the X-ray emission, and are 
consistent with the view that the apparent dichotomy be- 
tween the two BL Lac classes is mainly due to different cutoff 
energies. They would also bo consistent with the predictions 
of an inhomogoneous SSC model in which the main differ- 
ence is the orientation of the beaming axis, but this model 
appear not to be able to reproduce the observed large dif- 
ferences in cutoff energies. 



ACKNOWLEDGMENTS 

We thank Fabrizio Fiore, Gabriele Ghisellini, and Rita Sam- 
bruna for useful discussions, and Eric Perlman and Meg 
Urry for providing us with results in advance of publica- 
tion. This research has made use of the BROWSE pro- 
gram developed by the ESA/EXOSAT Observatory and by 
NASA/HEASARC. Wc acknowledge partial support from 
ESA/ESRIN, Frascati, where part of this work was done, 
while P.G. was the ESIS senior scientist and RR was a vis- 
iting scientist. 



Maraschi L., Fossati G., Tagliaferri G., Treves A., 1995, ApJ, 443, 
578 

Maraschi L., Ghisellini G., Tanzi E., Treves A., 1986, ApJ, 310, 
325 

Padovani P., Giommi P., 1995a, ApJ, 444, 567 
Padovani P., Giommi P., 1995b, MNRAS, 277, 1477 
Padovani P., Urry C. M., 1990, ApJ, 356, 75 
Perlman E. S. et al., 1995a, ApJS, in press 
Perlman E. S., Stocke J. T., 1993, ApJ, 406, 430 
Perlman E. S., Stocke J. T., Wang Q. D., Morris S. L., 1995b, 
ApJ, in press 

Rybicki G., Lightman A., 1979, Radiation Processes in Astro- 
physics. Wiley, New York 

Sambruna R. M., 1994, Ph. D. Thesis, SISSA 

Sambruna R. M., Barr P., Giommi P., Maraschi L., Tagliaferri 
G., Treves A., 1994, ApJ, 434, 468 

Stark A. A., Gammie C. F., Wilson R. W., Bally J., Linke R. A., 
Heiles C., Hurwitz M., 1992, ApJS, 79, 77 

Stickel M., Padovani P., Urry C. M., Fried J. W., Kuhr H., 1991, 
ApJ, 374, 431 

Stocke J. T., Licbcrt J., Schmidt G., Gioia I. M., Maccacaro T., 
Schild R. E., Maccagni D., Arp H. C., 1985, ApJ, 298, 619 

Urry C. M., Padovani P., 1995, PASP, 107, 803 

Urry C. M., Sambruna R. M., Worrall D. M., KoUgaard R. I., 
Feigelson E., Perlman E. S., Stocke J. T., 1996, ApJ, in press 

Veron-Cetty M.-P., Veron P., 1993, A Catalogue of Quasars and 
Active Nuclei {%^^ ed.), ESO Scientific Report n. 13 

White N. E., Giommi P., Angelini L., 1994, lAUC 6100 

NOTE ADDED IN PROOF 

After this paper was accepted for publication, it came to our 
knowledge that Lamer, Brunner & Staubert (1996), A&A, 
in press, have independently analysed most of the data pre- 
sented in this paper, reaching similar conclusions. 



REFERENCES 

Bregman J. N., 1990, A&A Review, 2, 125 

Brinkmann W., Siebert J., Reich W., Fiirst E., Rcicli P., Voges 
W., Trumper J., Wielebinski R., 1995, A&AS, 109, 147 

Brunner H., Lamer G., Worrall D. M., Staubert R., 1994, A&A, 
287, 436 

Ciliegi P., Bassani L., CaroU E., 1993, ApJS, 85, 111 

Ciliegi P., Bassani L., Caroli E., 1995, ApJ, 439, 80 

Comastri A., Molendi S., Ghisellini G., 1995, MNRAS, 277, 297 

Falomo R., Scarpa R., BersancUi M., 1994, ApJS, 93, 125 

Gear W. K., 1993, MNRAS, 264, L21 

Gear W. K. et al., 1986, ApJ, 304, 295 

George L M., Warwick R. S., Bromage G. E., 1988, MNRAS, 232, 
793 

Ghisellini G., Maraschi L., 1989, ApJ, 340, 181 

Giommi P., Ansari S. G., Micol A., 1995a, A&AS, 109, 267 

Giommi P., Barr P., Garilli B., Maccagni D., Pollock A. M. T., 

1990, ApJ, 356, 432 
Giommi P., Padovani P., White N. E., Angelini L., 1995b, in 

preparation 

Giommi P., Padovani, P., 1994, MNRAS, 268, L51 

Green R. F., Schmidt M., Liebert J., 1986, ApJS, 61, 305 

Hewitt A., Burbidge G., 1993, ApJS, 87, 451 

Jannuzi B. T., Green R. F., French H. 1993, ApJ, 404, 100 

Jannuzi B. T., Smith P. S., Elston R., 1994, ApJ, 428, 130 

KoUgaard R. I., 1994, Vistas in Astronomy, 38, 29 

Laurent-Muehleisen S. A., KoUgaard R. I., Moellenbrock G. A., 

Feigelson E. D., 1993, AJ, 106, 875 
Ledden J. E., O'DeU S. L., 1985, ApJ, 298, 630 



